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Forward 



The purpose of this book is to postulate some theories and test them numerically. 
Estimation is often a difficult task and it has wide application in social sciences and 
financial market. In order to obtain the optimum efficiency for some classes of 
estimators, we have devoted this book into three specialized sections: 



Part 1 . In this section we have studied a class of shrinkage estimators for shape 
parameter beta in failure censored samples from two-parameter Weibull distribution 
when some ’apriori’ or guessed interval containing the parameter beta is available in 
addition to sample information and analyses their properties. Some estimators are 
generated from the proposed class and compared with the minimum mean squared error 
(MMSE) estimator. Numerical computations in tenns of percent relative efficiency and 
absolute relative bias indicate that certain of these estimators substantially improve the 
MMSE estimator in some guessed interval of the parameter space of beta, especially for 
censored samples with small sizes. Subsequently, a modified class of shrinkage 
estimators is proposed with its properties. 



Part2. In this section we have analyzed the two classes of estimators for population 
median M Y of the study character Y using information on two auxiliary characters X and 
Z in double sampling. In this section we have shown that the suggested classes of 
estimators are more efficient than the one suggested by Singh et al (2001). Estimators 
based on estimated optimum values have been also considered with their properties. The 
optimum values of the first phase and second phase sample sizes are also obtained for the 
fixed cost of survey. 



Part3. In this section, we have investigated the impact of measurement errors on a family 
of estimators of population mean using multiauxiliary information. This error 
minimization is vital in financial modeling whereby the objective function lies upon 
minimizing over-shooting and undershooting. 



This book has been designed for graduate students and researchers who are active in the 
area of estimation and data sampling applied in financial survey modeling and applied 
statistics. In our future research, we will address the computational aspects of the 
algorithms developed in this book. 



The Authors 
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Abstract 

This paper is speculated to propose a class of shrinkage estimators for shape parameter 
ft in failure censored samples from two-parameter Weibull distribution when some ‘apriori ’ or 
guessed internal containing the parameter P is available in addition to sample information and 
analyses their properties. Some estimators are generated from the proposed class and compared 
with the minimum mean squared error (MMSE) estimator. Numerical computations in terms of 
percent relative efficiency and absolute relative bias indicate that certain of these estimators 
substantially improve the MMSE estimator in some guessed internal of the parameter space of p , 
especially for censored samples with small sizes. Subsequently, a modified class of shrinkage 
estimators is proposed with its properties. 



Key Words & Phrases: 

Two-parameter Weibull distribution, Shape parameter, Guessed internal, Shrinkage 
estimation technique, Absolute relative bias, Relative mean square error, Percent relative 
efficiency. 

2000 MSC: 62E17 

1. INTRODUCTION 

Identical rudiments subjected to identical environmental conditions will fail at different and 
unpredictable times. The ‘time of failure’ or ‘life length’ of a component, measured from some specified 
time until it fails, is represented by the continuous random variable X. One distribution that has been used 
extensively in recent years to deal with such problems of reliability and life-testing is the Weibull 
distribution introduced by Weibull(1939), who proposed it in connection with his studies on strength of 
material. 

The Weibull distribution includes the exponential and the Rayleigh distributions as special cases. 
The use of the distribution in reliability and quality control work was advocated by many authors following 
Weibull(1951), Lieblin and Zelen(1956), Kao(1958,1959), Berrettoni(1964) and Mann(1968 A). 
Weibull(1951) showed that the distribution is useful in describing the ‘wear-out’ or fatigue failures. 
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Kao(1959) used it as a model for vacuum tube failures while Lieblin and Zelen(1956) used it as a model for 
ball bearing failures. Mann(1968 A) gives a variety of situations in which the distribution is used for other 
types of failure data. The distribution often becomes suitable where the conditions for “strict randomness” 
of the exponential distribution are not satisfied with the shape parameter [S having a characteristic or 
predictable value depending upon the fundamental nature of the problem being considered. 

1.1 The Model 

Let Xi, x 2 , ..., x n be a random sample of size n from a two-parameter Weibull distribution, 
probability density function of which is given by : 

) = Pa p x p - X expj-(x/ a)^}; x > 0 ,a> 0,/? > 0 

(l.l) 

where a being the characteristic life acts as a scale parameter and ft is the shape parameter. 

The variable Y= In x follows an extreme value distribution, sometimes called the log-Weibull 
distribution [e.g. White(1969)], cumulative distribution function of which is given by : 

F(y) = 1-expj- exp^ 1 ^ — j j;-oo < y < oo,-oo < u < oo ,b > 0 

(1.2) 

where b = 1//? and u = In a are respectively the scale and location parameters. 

The inferential procedures of the above model are quite complex. Mann(1967 A,B, 1968 B) 
suggested the generalised least squares estimator using the variances and covariances of the ordered 
observations for which tables are available up to n = 25 only. 



1.2 Classical Estimators 

Suppose X\, Xn, ..., x m be the m smallest ordered observations in a sample of size n from Weibull 
distribution. Bain(1972) defined an unbiased estimator for b as 



(1.3) 
where 

(1.4) 
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y t ~ u 

and Vj = — ~ — are ordered variables from the extreme value distribution with u = 0 and b = 

A 

l.The estimator b u is found to have high relative efficiency for heavily censored cases. Contrary to this, 



the asymptotic relative efficiency of b u is zero for complete samples. 

Engelhardt and Bain(1973) suggested a general form of the estimator as 






where g is a constant to be chosen so that the variance of b g is least and K (g m n) is an unbiasing constant. 



U U g 

The statistic has been shown to follow approximately 7 - distribution with h degrees of freedom, 

b 



; 7 = 1,2 



where h = 2/ Var[ b g b . Therefore, we have 



R ip 1 _ _l f 2 1 ip r [( /?/2 ) + i> 

J p jp \h — 2J T{h!2) 



where [3 — — ~ — is an unbiased estimator of /3 with Var^pj= — — — and t = hb g having density 



m = 



i rp 



r(A/2)U 



£ ex p| -F f ( * /2H ; t > o 



The MMSE estimator of /?, among the class of estimators of the form C /3 ; C being a constant for 

A 

which the mean square error (MSE) of C /3 is minimum, is 

a h-4 

Pm = I ’ 



having absolute relative bias and relative mean squared error as 

arb {M = /F' 



7 




respectively. 



1.3 Shrinkage Technique of Estimation 

Considerable amount of work dealing with shrinkage estimation methods for the parameters of the 
Weibull distribution has been done since 1970. An experimenter involved in life-testing experiments 
becomes quite familiar with failure data and hence may often develop knowledge about some parameters of 
the distribution. In the case of Weibull distribution, for example, knowledge on the shape parameter ft can 
be utilised to develop improved inference for the other parameters. Thompson(1968 A,B) considered the 

problem of shrinking an unbiased estimator E, of the parameter C, either towards a natural origin E, or 



towards an interval 



|^,^jand suggested the shrunken estimators hE, + (1 — h)E,^ 



and 



/ 

V 



\ + £, ' 

~ 1 ~2 



, where 0 < h < 1 is a constant. The relevance of such type of shrunken 



estimators lies in the fact that, though perhaps they are biased, has smaller MSE than C, for C, in some 



interval around E or 

J A 



A ,+0 



, as the case may be. This type of shrinkage estimation of the Weibull 



parameters has been discussed by various authors, including Singh and Bhatkulikar(1978), Pandey(1983), 
Pandey and Upadhyay(1985,1986) and Singh and Shukla(2000). For example, Singh and 
Bhatkulikar(1978) suggested performing a significance test of the validity of the prior value of ft (which 
they took as 1). Pandey(1983) also suggested a similar preliminary test shrunken estimator for ft. 



In the present investigation, it is desired to estimate /? in the presence of a prior information 

available in the form of an interval (//, , // , ) and the sample information contained in ft . Consequently, 
this article is an attempt in the direction of obtaining an efficient class of shrunken estimators for the scale 
parameter f 3 . The properties of the suggested class of estimators are also discussed theoretically and 
empirically. The proposed class of shrunken estimators is furthermore modified with its properties. 

2. THE PROPOSED CLASS OF SHRINKAGE ESTIMATORS 

* 

Consider a class of estimators ft (p for ft in model (1.1) defined by 




p 



A + A) 




( \ 




q + w 


P 1 +P 2 




l 2 J 


{ J 





(P.q) 

( 2 . 1 ) 

where p and q are real numbers such that A 0 and q > 0, w is a stochastic variable which may in 

* 

particular be a scalar, to be chosen such that MSE of ft is minimum. 

* 

Assuming w a scalar and using result (1.6), the MSE of /3 is given by 

A 2 ) 2p r[(h/2) + 2p\ 



MSE ( 3 ( ^ =P 



{qA-l} 2 +w 2 A 2(p+1) 



h- 2. 



T{h/2) 



+ {qA-l}wA {p+1) 
(2.2) 



f 2 l 


^ r[(A/2) + 2/?] 


U-2j 


T{h!2) 



where 



A = 



P,+P; 

2p 



Minimising (2.2) with respect to w and replacing /3 by its unbiased estimator (3 , we get 



Pl + P2\_p[p P 



'1 ' H2 

2 



w = - 



P\ + Pi 

2 j 



\(p+ 1) 



-w(p). 



(2.3) 



where 



w(p) = 



h-2Y r[(h/2) + p\ 



2 ) T[{h!2) + 2 p\ 



(2.4) 



lies between 0 and 1, {i.e., 0 < w(p) < 1} provided gamma functions exist, i.e., p > {—hi 2) . 
Substituting (2.3) in (2.1) yields a class of shrinkage estimators for /? in a more feasible form as 



P 

(2.5) 

2.1 Non-negativity 



(A?) 



f h- 2^ 



V I J 



Mp) + q\ I I 1 - w (p)} ■ 



’ 1 1 P2 
2 
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Clearly, the proposed class of estimators (2.5) is the convex combination of {(/? — 2 )/ 1\ and 

(g(P[ +P 2 )/ 2} and hence P [pq) is always positive as {{h — '1)1 1] > 0 and q> 0. 

2.2 Unbiasedness 

If w(p) = 1, the proposed class of shrinkage estimators p ( ( turns into the unbiased estimator P , 
otherwise it is biased with 

Bias jPo>, ? ) | = P {^ A -!}[!- W (P)\ (2-6) 

and thus the absolute relative bias is given by 

ARB jP( M ) | = | -!}[!- HP)\ | • 

(2.7) 

The condition for unbiasedness that w(p) = 1, holds iff, censored sample size m is indefinitely 
large, i.e., m — > 00 . Moreover, if the proposed class of estimators $ (pqj turns into P then this case does not 
deal with the use of prior information. 

A more realistic condition for unbiasedness without damaging the basic structure of P fp q) and 
utilises prior information intelligibly can be obtained by (2.7). The ARB of P fp q) is zero when q = A 1 (or 

A = q~ [ ). 

2.3 Relative Mean Squared Error 

The MSE of the suggested class of shrinkage estimators is derived as 

MSE {p(^)} = P 2 {<?A- 1} 2 {1 -^)} 2 + 2 |j ( / 4 ) | ’ ( 2 - 8 ) 

and relative mean square error is therefore given by 

RMSE {P(p, ? , } = W A - 1} 2 {1 - w{p)Y + • 

(2.9) 

It is obvious from (2.9) that RMSE |j3 (/; q) j is minimum when q = A 1 (or A = q 1 ). 

2.4 Selection of the Scalar ‘p’ 
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The convex nature of the proposed statistic and the condition that gamma functions contained in 
w(p) exist, provides the criterion of choosing the scalar p. Therefore, the acceptable range of value of p is 
given by 

{p | 0<w(p)<\ and p>(-h/ 2)}, V n, m. (2.10) 



2.5 Selection of the Scalar ‘ q ’ 

It is pointed out that at q = A 1 , the proposed class of estimators is not only unbiased but renders 
maximum gain in efficiency, which is a remarkable property of the proposed class of estimators. Thus 
obtaining significant gain in efficiency as well as proportionately small magnitude of bias for fixed A or 
for fixed (Pj/P) and (P 2 /P), one should choose q in the vicinity of q = A 1 . It is interesting to note 

that if one selects smaller values of q then higher values of A leads to a large gain in efficiency (along 
with appreciable smaller magnitude of bias) and vice-versa. This implies that for smaller values of q. the 
proposed class of estimators allows to choose the guessed interval much wider, i.e., even if the 
experimenter is less experienced the risk of estimation using the proposed class of estimators is not higher. 
This is legitimate for all values of p. 



2.3 Estimation of Average Departure: A Practical Way of selecting q 

The quantity A = {(Pj + P 2 )/2p) , represents the average 

departure of natural origins fi x and (3 2 from the true value 
/? . But in practical situations it is hardly possible to get 
an idea about A . Consequently, an unbiased estimator of 
A is proposed, namely 

A p (p, +p 2 ) l r(/z/ 2 ) 

| 4 J T[(h!2) + \ ' 



( 2 . 12 ) 



In section 2.5 it is investigated that, if q= A 1 , the 
suggested class of estimators yields favourable results. 
Keeping in view of this concept, one may select q as 



q = A 



-i 



f 4 | r[(/z / 2) + 1 

IhPi+PdJ r(/i/2) 



(2.13) 



Here this is fit for being quoted that this is the 
criterion of selecting q numerically and one should 
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carefully notice that this doesn't mean q is replaced by 
(2.13) in fL., . 



3 . COMPARISION OF ESTIMATORS AND EMPIRICAL STUDY 

James and Stein(1961) reported that minimum MSE is a highly desirable property and it is 
therefore used as a criterion to compare different estimators with each other. The condition under which the 
proposed class of estimators is more efficient than the MMSE estimator is given below. 



MSE \J3 | does not exceed the MSE of MMSE estimator J3 M if - 



(i-Vg)^ 1 < a <(i +V g)^ 1 



where 



{l-w(p)} 2 [(h-2) (h- 4) 



Besides minimum MSE criterion, minimum bias is also important and therefore should be 



incorporated under study. Thus, ARB ' 



is less than ARB 



{Aujif- 



2 2 
1 ? x>cf l < A < < 1 h t \>q~ 



3.1 The Best Range of Dominance of A 

The intersection of the ranges of A in (3.1) and (3.2) gives the best range of dominance of A 
denoted by A Best . In this range, the proposed class of estimators is not only less biased than the MMSE 
estimator but is more efficient than that. The four possible cases in this regard are: 

(i) if ] 1 f i 1 < (l — *J(3 ) and i 1 H f n [•<(]+ VG ) then 

1 (h-2)[\-w(pj\) V ’ \ (A-2)[l-M<p)]J ' ’ 



’{ 1+ ( /i-2)[i-«<p)] r 



(ii) if 1- 



(h-2)[l-w(p)\\ 



1-VG and 1 + VG <0 + 



(h-2)\[-w(p)\ 



A B . is the same as defined in (3.1). 
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(iii) 



(iv) 



4 



l-VG < < 1 - 



(h ~ 2) [l - w(p)\ 



| and (l + VG j 



< 1 + 



{h-2)\Y-w(p)]\ 



then 



^ Best 



1 - 



(h-2)[l-w(pj\ 



\ c i 



if (i-Vg 



< 1 - 



(h-2)[l-w(p)\ 



\ i + Vg k 1 



| and |l + 



(h-2)[\-w(p)]\ 



< 



(i + Vg 



then 



A Best is the same as defined in (3.2). 



3.2 Percent Relative Efficiency 

A 

To elucidate the performance of the proposed class of estimators P {pq) w *th the MMSE 

A A A 

estimator P v/ , the Percent Relative Efficiencies (PREs) of P ( > with respect to P M have been computed 



by the formula: 



PRE 






2(h - 4) 



(h - 2)[(t)A - 1) 2 {l - w(p)} 2 (h - 4) + 2 {w(p)} 2 



x 100 



(3.5) The PREs of P ( ;) (/) with respect to p and ARBs of p (/) ^ for fixed n = 20 and different values 

of p,q,m Aj(= Pj/P) and A 2 (=p 2 /p) or A are compiled in Table 3.1 with corresponding values of h 
[which can be had from Engelhardt(1975)] and w(p). The first column in every m corresponds to PREs and 

A 

the second one corresponds to ARBs of P The last two rows of each set of q includes the range of 

dominance of A and A Best . The ARBs of P has also been given at the end of each set of table. 



M 
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Table 3.1 



A A A 

PREs of proposed estimator P with respect to MMSE estimator [) m and ARBs of P 



qi 


P = - 2 


Ail 


A2'l' 


/MS 


6 


8 


10 


12 




/is 


10.8519 


15.6740 


20.8442 


26.4026 








AnI- w{p) s 


0.1750 


0.3970 


0.5369 


0.6305 




0.1 


0.2 


0.15 


35.33 


0.7941 


40.20 




45.57 


0.4457 


50.60 






0.4 


0.6 


0.50 


42.62 


0.7219 


47.90 


0.5276 


53.49 




58.53 






0.4 


1.6 


1.00 


57.66 


0.6188 


63.18 


0.4522 


68.54 


0.3473 


72.99 






1.0 


2.0 


1.50 


82.21 


0.5156 


86.53 


0.3769 


89.95 


0.2894 


92.27 




0.25 


1.6 


2.4 


2.00 


126.15 


0.4125 


124.06 




120.83 


0.2315 


117.72 






2.0 


3.0 


2.50 


215.89 




187.20 


0.2261 


164.84 


0.1737 


149.86 






2.5 


3.5 


3.00 


438.90 




294.12 




222.82 


0.1158 


186.17 






3.5 


3.5 


3.50 


1154.45 




447.47 




282.42 




217.84 






3.8 


4.2 


4.00 


2528.52 




541.60 




310.07 




230.93 








Range of As 


(1.74, 




(1.70, 




(1.68, 




(1.66, 


(3.11, 










6.25) 




6.29) 




6.31) 


4.91) 


6.33) 


4.88) 




Aflest — > 












0.1 


0.2 


0.15 


38.21 


0.7632 


43.26 


0.5577 


48.75 


0.4284 


53.81 


0.3418 




0.4 


0.6 


0.50 


57.66 


0.6188 


63.18 


0.4522 


68.54 


0.3473 


72.99 


0.2771 




0.4 


1.6 


1.00 


126.15 


0.4125 


124.06 


0.3015 


120.83 


0.2315 


117.72 


0.1847 




1.0 


2.0 


1.50 


438.90 


0.2063 


294.12 


0.1507 


222.82 


0.1158 


186.17 




0.50 


1.6 


2.4 


2.00 


2528.52 




541.60 




310.07 




230.93 






2.0 


3.0 


2.50 


438.90 


IBS 


294.12 




222.82 


0.1158 


186.17 






2.5 


3.5 


3.00 


126.15 


0.4125 


124.06 




120.83 


0.2315 


117.72 


0.1847 




3.5 


3.5 


3.50 


57.66 


0.6188 


63.18 




68.54 


0.3473 


72.99 


0.2771 




3.8 


4.2 


4.00 


32.76 




37.45 




42.68 


0.4631 


47.65 


0.3695 






Range of As 


(0.87, 




(0.85, 




(0.84, 




(0.83, 












3.13) 




3.15) 




3.16) 




3.17) 






Aflest — > 












0.1 


0.2 


0.15 


41.45 


0.7322 


46.67 


0.5351 


52.25 






0.3279 




0.4 


0.6 


0.50 


82.21 


0.5156 


86.53 


0.3769 


89.95 


0.2894 


92.27 






0.4 


1.6 


1.00 






294.12 




222.82 


0.1158 


186.17 






1.0 


2.0 


1.50 


1154.45 


IB 


447.47 




282.42 




217.84 




0.75 


1.6 


2.4 


2.00 


126.15 


0.4125 


124.06 




120.83 


0.2315 


117.72 


0.1847 




2.0 


3.0 


2.50 


42.62 


0.7219 


47.90 


0.5276 


53.49 




58.53 


0.3233 




2.5 


3.5 


3.00 






24.58 


0.7537 


28.74 


0.5789 


32.94 


0.4619 




3.5 


3.5 


3.50 


12.51 




14.82 


0.9798 


17.67 


0.7525 








3.8 


4.2 


4.00 


8.27 




9.87 




11.90 


0.9262 










Range of As 






(0.57, 




(0.56, 


















nn 


2.10) 




2.11) 






mm m 




Aflest —> 










ARB of MMSE Estimators 




0.1463 


0.1061 





14 

















Table 3.1 continued ... 







A 2' l ' 


h-> 








Axt 




0.1 


0.2 


0.15 




0.4 


0.6 


0.50 




0.4 


1.6 


1.00 




1.0 


2.0 


1.50 


0.25 


1.6 


2.4 


2.00 




2.0 


3.0 


2.50 




2.5 


3.5 


3.00 




3.5 


3.5 


3.50 




3.8 


4.2 


4.00 




Range of A^ 




Aflest 




0.1 


0.2 


0.15 1 




0.4 


0.6 


0.50 




0.4 


1.6 


1.00 




1.0 


2.0 


1.50 


0.50 


1.6 


2.4 


2.00 




2.0 


3.0 


2.50 




2.5 


3.5 


3.00 




3.5 


3.5 


3.50 




3.8 


4.2 


4.00 




Range of A— » 




Aflest — > 




0.1 


0.2 


0.15 




0.4 


0.6 


0.50 




0.4 


1.6 


1.00 




1.0 


2.0 


1.50 


0.75 


1.6 


2.4 


2.00 




2.0 


3.0 


2.50 




2.5 


3.5 


3.00 




3.5 


3.5 


3.50 




3.8 


4.2 


4.00 




Range of A^ 




Aflest — > 


ARB of MMSE Estimator — > 



6 

10.8519 

0.7739 



P =-l 
8 

15.6740 

0.8537 



10 

20.8442 

0.8939 



12 

26.4026 

0.9180 




03.38 
10.98 
20.43 
126.91 0.0565 

129.23 0.0000 

126.91 0.0565 

120.43 0.1130 

110.98 0.1696 

100.00 0.2261 
(0.00, (0.00, 
4.00) 4.00) 

(0.00, 4.00) 



105.05 

115.99 

126.91 

128.65 

120.43 

105.60 

88.71 

72.93 

59.57 

( 0 . 00 , 

2.67) 

(0.00 



0.2259 








103.21 

109.79 

115.89 

116.82 

112.32 

103.55 

92.40 

80.65 

69.50 

( 0 . 00 , 

2.67) 

( 0 . 00 . 



0.1463 





102.31 0.0942 

106.91 0.0663 

111.05 0.0265 

111.67 0.0133 

108.65 0.0531 

102.55 0.0929 

94.37 0.1327 

85.17 0.1725 

75.85 0.2123 

(0.00, (0.00, 

2.67) 2.67) 

(0.00, 2.67) 



0.1061 



103.73 0.0615 

105.27 0.0513 

106.56 0.0410 

107.59 0.0308 

108.34 0.0205 

108.79 0.0103 

108.94 0.0000 

(0.00, (0.00, 
8.00) 8.00) 
(0.00, 8.00) 



1.20 
3.73 

106.56 0.0410 

108.34 0.0205 

108.94 0.0000 

108.34 0.0205 

106.56 0.0410 

103.73 0.0615 

100.00 0.0820 
(0.00, (0.00, 
4.00) 4.00) 

(0.00, 4.00) 



101.77 

105.27 

108.34 

108.79 

106.56 

101.96 

95.59 

88.13 

80.24 

( 0 . 00 , 

2.67) 

( 0 . 00 . 



0.0820 
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Table 3.1 continued ... 



qi 


P =1 




A 2' l ' 


m — > 


6 


8 


10 


12 




/*-> 


10.8519 


15.6740 


20.8442 


26.4026 








Al >v(/7) — > 


0.6888 


0.7737 


0.8251 


0.8779 




0.1 


0.2 


0.15 


99.00 




97.51 


0.2178 


97.21 




mm 


0.1175 




0.4 


0.6 


0.50 


106.26 




103.17 




101.80 




mm 






0.4 


1.6 


1.00 


117.09 




111.34 


0.1697 


108.25 




11 


0.0916 




1.0 


2.0 


1.50 


128.15 




119.34 


0.1415 


114.39 


0.1093 


mm 


0.0763 


0.25 


1.6 


2.4 


2.00 


138.88 




126.79 


0.1132 


119.95 


0.0875 


mm 


0.0611 




2.0 


3.0 


2.50 


148.56 




133.27 




124.67 


0.0656 


BB 


0.0458 




2.5 


3.5 


3.00 


156.33 




138.31 




128.27 


0.0437 


117.53 






3.5 


3.5 


3.50 


161.41 




141.52 




130.54 


0.0219 


118.72 


0.0153 




3.8 


4.2 


4.00 


163.17 




142.63 




131.31 


0.0000 


119.12 








Range of A — » 






(0.30, 




(0.36, 


(0.00, 
















M™ 


7.70) 




7.64) 


8.00) 


7.76) 












(0.36, 7.64) 






0.1 


0.2 


0.15 


102.07 


0.2879 


99.92 




99.18 


0.1618 


100.49 


0.1130 




0.4 


0.6 


0.50 


117.09 


0.2334 


111.34 


0.1697 


108.25 


0.1312 


106.18 


0.0916 




0.4 


1.6 


1.00 


138.88 


0.1556 


126.79 


0.1132 


119.95 


0.0875 


113.00 


0.0611 




1.0 


2.0 


1.50 


156.33 




138.31 




128.27 


0.0437 


117.53 


0.0305 


0.50 


1.6 


2.4 


2.00 


163.17 




142.63 




131.31 




119.12 






2.0 


3.0 


2.50 


156.33 




138.31 




128.27 


0.0437 


117.53 






2.5 


3.5 


3.00 


138.88 


0.1556 


126.79 


0.1132 


119.95 


0.0875 


113.00 






3.5 


3.5 


3.50 




0.2334 


111.34 


0.1697 




0.1312 


106.18 






3.8 


4.2 


4.00 




0.3112 


95.12 


0.2263 


95.25 


0.1749 


97.90 








Range of A — » 






































3.88) 


3.34) 




Aaest — > 












0.1 


0.2 


0.15 


105.20 


0.2762 


102.36 




101.15 


Mjggg 


101.75 






0.4 


0.6 


0.50 


128.15 


0.1945 


119.34 


0.1415 


114.39 


HE 


109.82 






0.4 


1.6 


1.00 


156.33 




138.31 




128.27 


0.0437 


117.53 






1.0 


2.0 


1.50 


161.41 




141.52 




130.54 


0.0219 


118.72 




0.75 


1.6 


2.4 


2.00 


138.88 


0.1556 


126.79 


0.1132 


119.95 


0.0875 


113.00 






2.0 


3.0 


2.50 




0.2723 


m 




101.80 


0.1531 


102.17 






2.5 


3.5 


3.00 


77.96 


0.3891 


an 


0.2829 


82.50 


0.2187 


88.98 






3.5 


3.5 


3.50 


57.31 




61.51 


0.3678 


65.66 


0.2843 


75.76 






3.8 


4.2 


4.00 


42.96 


0.6225 


47.58 


0.4526 


52.22 


0.3499 


63.80 








Range of A — » 










(0.12, 




(0.08, 




















2.55) 




2.59) 


2.23) 




Aflest — > 






(0.52,2.14) 




ARB of MMSE Estimator— » 


0.2259 


0.1463 


0.1061 


0.0820 
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Table 3.1 continued ... 



qi 


P =2 


Ail 


A2'i' 


in— > 


6 


8 


10 


12 




/*-> 


10.8519 


15.6740 


20.8442 


26.4026 








Al w ( p ) — > 


0.3131 


0.4385 


0.5392 


0.6816 




0.1 


0.2 


0.15 


48.51 


0.6612 


45.00 




45.90 






mgm 




0.4 


0.6 


0.50 


57.95 


0.6011 


53.31 


0.4913 


53.85 




68.81 


mm 




0.4 


1.6 


1.00 


76.84 


0.5152 


69.55 


0.4211 


68.94 


0.3456 




IB 




1.0 


2.0 


1.50 


106.11 


0.4293 


93.70 


0.3509 


90.35 




101.08 


0.1990 


0.25 


1.6 


2.4 


2.00 


154.14 


0.3435 


130.87 


0.2808 


121.15 




122.65 


0.1592 




2.0 


3.0 


2.50 


237.92 


0.2576 


189.27 


0.2106 


164.85 


0.1728 


147.06 


0.1194 




2.5 


3.5 


3.00 


388.87 


0.1717 


277.82 


0.1404 


222.08 


0.1152 


171.43 


0.0796 




3.5 


3.5 


3.50 


627.92 


0.0859 


386.26 




280.49 




190.36 


0.0398 




3.8 


4.2 


4.00 


789.74 


0.0000 






307.45 




197.63 


0.0000 






Range of A — » 


( 1 - 41 , 


( 2 . 68 , 






( 1 . 68 , 




( 1 . 47 , 


( 2 . 97 , 










6 . 59 ) 


5 . 32 ) 






6 . 32 ) 


4 . 92 ) 


6 . 53 ) 


5 . 03 ) 




Afiest — > 


( 2 . 68 , 5 . 32 ) 






( 2 . 97 , 5 . 03 ) 




0.1 


0.2 


0.15 


52.26 


0.6354 


48.32 


0.5194 


49.09 


0.4262 


63.91 






0.4 


0.6 


0.50 


76.84 


0.5152 


69.55 


0.4211 


68.94 


0.3456 


83.20 






0.4 


1.6 


1.00 


154.14 


0.3435 


130.87 




121.15 




122.65 






1.0 


2.0 


1.50 


388.87 


0.1717 


277.82 




222.08 


0.1152 


171.43 


0.0796 


0.50 


1.6 


2.4 


2.00 


789.74 




444.03 




307.45 




197.63 


0.0000 




2.0 


3.0 


2.50 


388.87 


0.1717 


277.82 




222.08 


0.1152 


171.43 


0.0796 




2.5 


3.5 


3.00 


154.14 


0.3435 


130.87 




121.15 




122.65 


0.1592 




3.5 


3.5 


3.50 


76.84 


0.5152 


69.55 


0.4211 


68.94 


0.3456 


83.20 


0.2388 




3.8 


4.2 


4.00 


45.14 


0.6869 


42.00 


0.5615 


42.99 




57.36 


0.3184 






Range of A — » 


( 0 . 71 , 




( 0 . 80 , 




( 0 . 84 , 




( 0 . 74 , 


( 1 . 49 , 










3 . 29 ) 




3 . 20 ) 




3 . 16 ) 




3 . 26 ) 


2 . 51 ) 




Aflest — > 








( 1 . 49 , 2 . 51 ) 




0.1 


0.2 


0.15 


56.45 




52.00 


0.4983 


52.60 




67.54 


0.2826 




0.4 


0.6 


0.50 


106.11 


0.4293 


93.70 




90.35 




101.08 






0.4 


1.6 


1.00 


388.87 


0.1717 


277.82 




222.08 


0.1152 


171.43 






1.0 


2.0 


1.50 


627.92 




386.26 




280.49 




190.36 




0.75 


1.6 


2.4 


2.00 


154.14 


0.3435 


130.87 




121.15 




122.65 


0.1592 




2.0 


3.0 


2.50 


57.95 




53.31 


0.4913 


53.85 




68.81 


0.2786 




2.5 


3.5 


3.00 


29.50 


0.8587 


27.83 




28.97 




41.00 






3.5 


3.5 


3.50 


17.73 


1.1163 


16.90 


0.9125 


17.83 


0.7488 


26.50 


0.5175 




3.8 


4.2 


4.00 


11.79 


1.3739 


11.30 




12.01 


0.9216 


18.33 


0.6369 






Range of A — » 


( 0 . 47 , 




( 0 . 53 , 




( 0 . 56 , 




( 0 . 49 , 


( 0 . 99 , 










2 . 20 ) 




2 . 13 ) 


1 . 68 ) 


2 . 11 ) 




2 . 18 ) 


1 . 68 ) 




Afiest — > 










ARB of MMSE Estimator— > 


0.2259 


0.1463 


0.1061 


0.0820 
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3000.00 
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n = 20 
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It has been observed from Table 3.1, that on keeping m, p, q fixed, the relative efficiencies of the 
proposed class of shrinkage estimators increases lip to A = q~ l , attains its maximum at this point and then 
decreases symmetrically in magnitude, as A increases in its range of dominance for all n, p and q. On the 
other hand, the ARBs of the proposed class of estimators decreases up to A = q~\ the estimator becomes 
unbiased at this point and then ARBs increases symmetrically in magnitude, as A increases in its range of 
dominance. Thus it is interesting to note that, at q = A~' , the proposed class of estimators is unbiased with 
largest efficiency and hence in the vicinity of q = A' 1 also, the proposed class not only renders the massive 
gain in efficiency but also it is marginally biased in comparison of MMSE estimator. This implies that q 
plays an important role in the proposed class of estimators. The following figure illustrates the discussion. 

Figure 3.1 

The effect of change in censored sample size m is also a matter of great interest. For fixed p. q and 
A , the gain in relative efficiency diminishes, and ARB also decreases, with increment in m. Moreover, it 
appears that to get better estimators in the class, the value of w(p) should be as small as possible in the 
interval (0,1]. Thus, to choose p one should not consider the smaller values of w(p) in isolation, but also the 
wider length of the interval of A. 
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4. MODIFIED CLASS OF SHRINKAGE ESTIMATORS AND ITS PROPERTIES 

The proposed class of estimators P ( ( is not uniformly better than p . It will be better if |3j and 

P 2 are in the vicinity of true value p . Thus, the centre of the guessed interval (p, + P 2 )/ 2 is of much 
importance in this case. If we partially violate this, i.e., only the centre of the guessed interval is not of 
much importance, but the end points of the interval p and p 2 are itself equally important then we can 
propose a new class of shrinkage estimators for the shape parameter P by using the suggested class 

P(m) as 

'P, , if t>[(h- 2)/pJ 

P ( P , q) = ' ^-1 w (p) + q[ — {l - HP)} , if [(h ~ 2)/P 2 ] < t < [(/? - 2)/pJ 

V t J V z J 

JS 2 , if <<[(«-2)/P 2 ] 

(4.1) 

which has 
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This modified class of shrinkage estimators is proposed in accordance with Rao(1973) and it 
seems to be more realistic than the previous one as it deals with the case where the whole interval is taken 
as apriori information. 



5. NUMERICAL ILLUSTRATIONS 

The percent relative efficiency of the proposed estimator p ( q) with respect to MMSE 



estimator P m has been defined as 



PRE © (M) ,P B = 



MSEp 



MSE p 






x 100 



(5.1) 



and it is obtained for n = 20 and different values of p, q, m, A l and A, (or A). The findings are 
summarised in Table 5.1 with corresponding values of h and w(p). 



Table 5.1 

^ A 

PREs of proposed estimator p ( q) with respect to MMSE estimatorp m 













n 


= 20 












qi 


P -» 




- 


1 




1 




A2>1' 


m — » 


6 


8 


10 


12 


6 


8 


10 


12 


h 


10.8519 


15.6740 


20.8442 


26.4026 


10.8519 


15.6740 


20.8442 


26.4026 


A-l w(p)^> 


0.7739 


0.8537 


0.8939 


0.9180 


0.6888 


0.7737 


0.8251 


0.8779 




EB 


0.3 


0.25 


mm 


41.39 


34.91 


30.59 


49.84 




34.66 


31.15 






0.6 


0.50 


MM 




67.45 


63.17 




79.57 


65.63 


61.55 


■ 


19 


0.9 


0.75 


261.72 


227.42 


203.08 


172.06 


227.59 


191.97 


172.31 


156.69 


0.25 


0.8 


1.2 






426.98 


342.54 


286.06 


454.93 


355.31 


293.42 


262.79 




ID 


1.5 


1.25 


649.95 




375.91 


314.98 


636.21 




427.74 


353.74 






1.8 


1.50 


268.31 


189.82 




125.21 






168.38 






IS 


2.0 


1.75 




53.66 




31.38 


82.35 






31.74 






0.3 


0.25 


50.84 


41.32 


34.76 


30.39 


49.90 


40.03 


34.45 








0.6 


0.50 


120.81 


82.01 


67.97 


63.49 


118.31 


81.13 


66.48 




■ 


is 


0.9 


0.75 


298.17 


253.12 


221.74 


184.38 


271.73 


225.47 




173.57 


0.50 


0.8 


1.2 




642.86 


473.19 


368.65 


303.15 


583.65 


433.16 




292.64 




ID 


1.5 


1.25 


626.09 


435.87 


345.16 


289.53 


658.77 


481.87 




317.87 




H 


1.8 


1.50 




175.97 




118.43 


264.16 


191.09 


152.66 


124.73 




IB 


2.0 


1.75 


78.41 


52.66 


39.39 


31.11 


79.96 


53.72 




31.36 






0.3 


0.25 


50.89 


41.24 


34.60 


30.19 


49.97 


39.95 


34.23 








0.6 


0.50 


124.02 


83.01 


68.50 


63.81 


122.74 


82.68 


67.32 


62.50 


■ 




0.9 


0.75 


339.92 


282.24 


242.46 


197.73 


325.66 


266.36 


229.58 


192.68 


0.75 


0.8 


1.2 


1.00 


723.50 


510.42 


389.34 


316.87 


710.96 


504.67 


388.35 


317.53 




1.0 


1.5 


1.25 


566.19 


392.47 


312.16 


263.77 


597.64 


421.61 


337.17 


278.26 
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1.2 


1.8 


1.50 


224.67 


161.95 


131.14 


111.81 


233.41 


169.19 


136.65 


114.63 


1.5 


2.0 


1.75 


76.05 


51.59 


38.85 


30.83 


76.93 


52.14 


39.17 


30.95 
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Table 5.1 continued ... 





P -> 




- 


2 




2 




A,i 


A2>1' 


m — » 


6 


8 


10 


12 


6 


8 


10 


12 




It 


10.8519 


15.6740 


20.8442 


26.4026 


10.8519 


15.6740 


20.8442 


26.4026 








A-l w(p)— > 


0.7739 


0.8537 


0.8939 


0.9180 


0.6888 


0.7737 


0.8251 


0.8779 






0.3 




46.04 


34.18 


30.92 


30.53 


46.77 


34.81 


30.96 


31.23 






0.6 




92.48 


72.59 


59.44 


53.42 


98.00 


73.36 


59.48 


54.88 


1 




0.9 


| 


106.83 


95.44 


92.75 


90.11 


128.68 


102.24 


93.16 


100.45 


0.25 


0.8 


1.2 


■SB 


145.02 


131.16 


126.15 


122.15 


191.47 


145.23 


126.97 


144.22 




ID 


1.5 


1.25 


220.29 


243.10 


282.54 


320.74 


305.32 


273.81 


284.60 


368.42 




IE 


1.8 


1.50 


208.14 


211.32 


202.36 


179.81 


250.20 


220.57 


202.56 


175.49 




m 


2.0 


1.75 










84.21 












0.3 


MEM 


46.28 


34.31 






46.95 


34.91 










0.6 




103.18 


76.82 


61.54 






77.31 


61.57 




M 




0.9 


| 


157.81 


135.64 




118.59 




142.94 


127.44 


128.23 


0.50 


0.8 


1.2 


mem 


267.16 








331.58 


246.71 








ID 


1.5 


1.25 


445.44 




448.55 


438.38 




467.49 


449.42 


432.21 




IE 


1.8 


1.50 






198.56 


163.98 


298.93 


238.16 








m 


2.0 


1.75 


84.92 


57.28 


42.13 


32.67 


84.44 




42.12 


32.44 






0.3 


MEM 


46.50 


34.43 


30.78 


29.92 


47.13 


34.99 


30.82 


30.50 






0.6 




114.64 


81.04 


63.59 


56.13 


116.87 


81.23 


63.61 


57.24 


■ 




0.9 


EH 


247.11 


202.90 


181.31 


160.85 


266.60 


209.00 


181.65 


167.34 


0.75 


0.8 


1.2 


WEM 


543.26 


418.40 


345.15 


293.90 


596.79 


430.93 


345.67 


302.22 




ID 


1.5 


1.25 


704.42 


541.77 


447.06 


381.03 


696.36 


532.12 


446.25 


358.48 




IE 


1.8 


1.50 


280.39 


203.46 


160.74 


132.95 


269.47 


199.82 


160.55 


129.07 




m 


2.0 


1.75 


81.39 


54.49 


40.40 


31.66 


80.35 


54.26 


40.39 


31.52 



It has been observed from Table 5.1 that likewise (3 (/) q) the PRE of [], „ q) with respect to (3 m 
decreases as censoring fraction (min) increases. For fixed m, p and q the relative efficiency increases lip to 
a certain point of A , procures its maximum at this point and then starts decreasing as A increases. It 
seems from the expression in (4.3) that the point of maximum efficiency may be a point where either any 
one of the following holds or any two of the following holds or all the following three holds- 

(i) the lower end point of the guessed interval, i.e., (3j coincides exactly with the true value (3 , i.e., 
A^l. 

(ii) the upper end point of the guessed interval, i.e., P 2 departs exactly two times from the true value 
p, i.e., A 2 = 2. 

(iii) A = q 1 

This leads to say that on contrary to P (p? ), there is much importance of A, and A 2 in addition to A. 
The discussion is also supported by the illustrations in Table 5.1. As well, the range of dominance of 
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average departure A is smaller than that is obtained for p {/j but this does not humiliate the merit of 
P ( q) because still the range of dominance of A is enough wider. 

6. CONCLUSION AND RECOMMENDATIONS 

It has been seen that the suggested classes of shrunken estimators have considerable gain in 
efficiency for a number of choices of scalars comprehend in it, particularly for heavily censored samples, 
i.e., for small m. Even for buoyantly censored samples, i.e., for large m, so far as the proper selection of 
scalars is concerned, some of the estimators from the suggested classes of shrinkage estimators are more 
efficient than the MMSE estimators subject to certain conditions. Accordingly, even if the experimenter has 
less confidence in the guessed interval (P j , p 2 ) of (3, the efficiency of the suggested classes of shrinkage 
estimators can be increased considerably by choosing the scalars p and q appropriately. 

While dealing with the suggested class of shrunken estimators p (/J q) it is recommended that one 
should not consider the substantial gain in efficiency in isolation, but also the wider range of dominance of 
A , because enough flexible range of dominance of A will leads to increase the possibility of getting 
better estimators from the proposed class. Thus it is recommended to use the proposed class of shrunken 
estimators in practice. 
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Abstract: 

In this paper we have suggested two classes of estimators for population median M Y of the study character 
Y using information on two auxiliary characters X and Z in double sampling. It has been shown that the 
suggested classes of estimators are more efficient than the one suggested by Singh et al (2001). Estimators 
based on estimated optimum values have been also considered with their properties. The optimum values 
of the first phase and second phase sample sizes are also obtained for the fixed cost of survey. 
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1. INTRODUCTION 

In survey sampling, statisticians often come across the study of variables which have highly skewed 
distributions, such as income, expenditure etc. In such situations, the estimation of median deserves special 
attention. Kuk and Mak (1989) are the first to introduce the estimation of population median of the study 
variate Y using auxiliary information in survey sampling. Francisco and Fuller (1991) have also 
considered the problem of estimation of the median as part of the estimation of a finite population 
distribution function. Later Singh et al (2001) have dealt extensively with the problem of estimation of 
median using auxiliary information on an auxiliary variate in two phase sampling. 

Consider a finite population U={l,2,...,i,...,N}. Let Y and X be the variable for study and auxiliary 
variable, taking values Y; and X; respectively for the i-th unit. When the two variables are strongly related 
but no information is available on the population median M x of X, we seek to estimate the population 
median M Y of Y from a sample S m , obtained through a two-phase selection. Permitting simple random 
sampling without replacement (SRSWOR) design in each phase, the two-phase sampling scheme will be as 
follows: 



(i) The first phase sample S n (S n cU) of fixed size n is drawn to observe only X in order to 
furnish an estimate of M x . 

(ii) Given S n , the second phase sample S m (S m <zS n ) of fixed size m is drawn to observe Y 
only. 

Assuming that the median M x of the variable X is known, Kuk and Mak (1989) suggested a ratio estimator 
for the population median M Y of Y as 
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(1.1) 



M, = My H- 
M 



x_ 

x 



where M Y and M x are the sample estimators of M Y and M x respectively based on a sample S m of size 
m. Suppose that y ( i,, y (2 ), . . y( m ) are the y values of sample units in ascending order. Further, let t be an 
integer such that Y (t) < M Y <Y (t -i , and let p=t/m be the proportion of Y, values in the sample that are less 
than or equal to the median value M Y , an unknown population parameter. If p is a predictor of p, the 

sample median M Y can be written in terms of quantities as Q y ( p ) where p = 0.5 . Kuk and Mak 
(1989) define a matrix of proportions (P;j(x,y)) as 





Y<M y 


Y>M y 


Total 


X < M x 


Pn(x,y) 


P 2 i(x,y) 


P i(x,y) 


X > M x 


Pn(x,y) 


P 22 (x,y) 


P 2 (x,y) 


Total 


P r(x,y) 


P 2 -(x,y) 


1 



and a position estimator of M y given by 



M r M =Qr(p r ) 



where 



Py 




' m x P\ i (x, v) | (m-m x )p n (x,y) 

v p.i(x,y) p,_(x,y) 

m.xPi i (*» y) + ( m ~ m ,-)Pii (*» y) ' 

m y 



A 



J 



(1.2) 



with Pj ( X , y) being the sample analogues of the Pi,(x,y) obtained from the population and m x the number 
of units in S m with X < M x . 



Let F ya ( v) and F YB ( v) denote the proportion of units in the sample S m with X < M x , and X>M x , 
respectively that have Y values less than or equal to y. Then for estimating M Y , Kuk and Mak (1989) 
suggested the 'stratification estimator' as 

M Y {St) = inf {y : F Y (y) > 0.5} (1.3) 



where F y ( v) = 



1 

2 



p (r) 

' YA 



+ F, 



YB 



(v) 



It is to be noted that the estimators defined in (1.1), (1.2) and (1.3) are based on prior knowledge of the 
median M x of the auxiliary character X. In many situations of practical importance the population median 
M x of X may not be known. This led Singh et al (2001) to discuss the problem of estimating the 
population median M Y in double sampling and suggested an analogous ratio estimator as 



= My 



M\ 



M 



X 



(1.4) 
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where M \ is sample median based on first phase sample S n . 

Sometimes even if M x is unknown, information on a second auxiliary variable Z, closely related to X but 
compared X remotely related to Y, is available on all units of the population. This type of situation has 
been briefly discussed by, among others, Chand (1975), Kiregyera (1980, 84), Srivenkataramana and Tracy 
(1989), Sahoo and Sahoo (1993) and Singh (1993). Let M z be the known population median of Z. 

Defining 



My 

M v 



-1 



= 



'Mj L _ 1 ^ 

y M X J 



, e 2 






Mz_ 

yM z 



\ 

1 and e 4 

J 



M 7 



1 



such that E(ek)s0 and I ek I <1 for k=0, 1 ,2,3; where M 2 and M\ are the sample median estimators based 
on second phase sample S m and first phase sample S n . Let us define the following two new matrices as 



and 





Z < M z 


Z > M z 


Total 


X < M x 


Pn(x,z) 


P 2 i(x,z) 


Pi(x,z) 


X > M x 


P !2 (x,z) 


P 22 (x,z) 


P 2 (x,z) 


Total 


P r (x,z) 


P 2 -(x,z) 


1 





Z < M z 


Z > M z 


Total 


Y<M y 


Pn(y,z) 


P 2 i(y,z) 


Pi(y,z) 


Y>M y 


Pi 2 (y,z) 


P 22 (y,z) 


P 2 (y,z) 


Total 


Pr(y,z) 


P 2 -(y,z) 


l 



Using results given in the Appendix-1, to the first order of approximation, we have 

E(e 0 2 ) = (4m)- 1 {M Y f Y (M Y )}- 2 , 

E(e > 2) = (^hP) (4m) 1 {M x fx(M x )}- 2 , 

E(e 2 2 ) = (^f)(4n)- 1 {M x f x (M x )}- 2 , 

E(e 3 2 ) = (^f) (4m)- 1 {M z f z (M z )}- 2 , 

E(e 4 2 ) = (^)(4n)- 1 {M z f z (M z )}- 2 , 

Efeoe,) = (4m)' 1 {4P u (x,y)-l} {M x M Y f x (M x )f Y (M Y )} ‘, 

E(e 0 e 2 ) = (Pf) (4n)- 1 {4P„(x,y)-l}{M x M Y f x (M x )f Y (M Y )}- 1 , 

E(e 0 e 3 ) = (^p) (4m)" 1 {4P i i(y,z)- 1 } {M Y M z f Y (M Y )f z (M z ) } 1 , 

E(e 0 e 4 ) = (Pp) (4n)- 1 {4P u (y,z)-l}{M Y M z f Y (M Y )f z (M z )}- 1 , 

E( ei e 2 ) = (Pp) (4n)-’{M x f x (M x )}- 2 , 

E(eie 3 ) = (ppp (4m) _l {4Pn(x,z)-l} (MxMzf^M^f^Mz)}" 1 , 
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E( ei e 4 ) = (4n)- 1 {4P 11 (x,z)-l}{M x M z fx(M x )fz(M z )}- 1 , 

E(e 2 e 3 ) = (4n)- 1 {4P u (x,z)-l}{M x M z fx(M x )f z (M z )}- 1 , 

E(e 2 e 4 ) = (4n)- 1 {4P 11 (x,z)-l}{M x M z fx(M x )f z (M z )}- 1 , 

E(e 3 e 4 ) = (^)(4n)- 1 (f z (M z )M z )- 2 

where it is assumed that as N— >co the distribution of the trivariate variable (X,Y,Z) approaches a continuous 
distribution with marginal densities f x (x), f Y (y) and f z (z) for X, Y and Z respectively. This assumption 
holds in particular under a superpopulation model framework, treating the values of (X, Y, Z) in the 
population as a realization of N independent observations from a continuous distribution. We also assume 
that f Y (M Y ), f x (M x ) and f z (M z ) are positive. 



Under these conditions, the sample median M Y is consistent and asymptotically normal (Gross, 1980) with 
mean M Y and variance 



' N — m ^ 

v N , 

In this paper we have suggested a class of estimators for M Y using information on two auxiliary variables X 
and Z in double sampling and analyzes its properties. 

2. SUGGESTED CLASS OF ESTIMATORS 




Motivated by Srivastava (1971), we suggest a class of estimators of M Y of Y as 

g = 



\m y {s) -.M y {s) =M Y g(u,v )} 



(2.1) 



where II = 



M, 



M 



,v = ■ 



M\ 

M, 



and g(u,v) is a function of u and v such that g(l,l)=l and such that it satisfies 



x z 

the following conditions. 



1. Whatever be the samples (S n and S m ) chosen, let (u,v) assume values in a closed convex sub- 
space, P, of the two dimensional real space containing the point (1,1). 

2. The function g(u,v) is continuous in P, such that g(l,l)=l. 

3. The first and second order partial derivatives of g(u,v) exist and are also continuous in P. 

Expanding g(u,v) about the point (1,1) in a second order Taylor's series and taking expectations, it is found 
that 



e[m y ^ 8) )=m y +00T 1 ) 

so the bias is of order n _1 . 

Using a first order Taylor's series expansion around the point (1,1) and noting that g(l,l)=l, we have 
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My (g) =My[l + e 0 +(e 1 -e 2 )g 1 (l,l) + e 4 g 2 (l,l) + o(H H )] 



or 



(m[ s) -M Y )=M Y [e Q +(e x - e 2 )g l (l,\) + e 4 g 2 (\,l)] (2.2) 



where gi(l,l) and g 2 ( 1,1) denote first order partial derivatives of g(u,v) with respect to u and v respectively 
around the point (1,1). 

Squaring both sides in (2.2) and then taking expectations, we get the variance of M Y ^ a ' > to the first degree 
of approximation, as 



Var\M 



(i My (g) ): 



1 



4(/r(M y )) 2 LU N) 



1 1 



+ 



V m n ) 



A + 



[^n N ) 



B 



(2.3) 



where 



A = 



B = 



M Y f Y {M Y ) ' 
M x f x {M x \ 



gi(U) 



M Y fy{M Y )\ ( ^ 

M z f z {M z )J Sz 1,1 



M Y fy(M Y ) 
M x f x [M x ) 

M Y f Y {M Y ) 
M zfz^M z) . 



g x (\,\)+2{4P n {x,y)-\) 



g 2 (l,l) + 2(4P 1I (j;,z)-l) 



(2.4) 



(2.5) 



/v / \ 

The variance of M Y ^ in (2.3) is minimized for 

'M x / x {M x y 



gi(U) = -| 
g 2 (U) = -| 



Myfy{M Y ) 

M zfziM z) 

Myfy{My J y 



(4Pn(x,v)-l) 



(4P n (y,z)-l) 



Thus the resulting (minimum) variance of M Y ^° ^ is given by 



(2.6) 



min. Var M 






4(/ y (M 7 )) 2 L^» N 



1 1 






i i 




llVT> z ) _ l) 
(2.7) 



Now, we proved the following theorem. 
Theorem 2.1 - Up to terms of order n. 



Var(M Y g )> 



1 



4 (f y {My)f 



in N 



T - I](4P, , (x,y)~ If - (i - £)(4P, , ( v, z) - 1) 2 



with equality holding if 
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&i(U) = - 

g 2 (U) = -| 



Mj x {M x y 

Myfy{M Y ) 

M z f z (M z ) 
Myf Y {M Y ) 



(4P n (x,y)-l) 



(4P n (y,z)-l) 



It is interesting to note that the lower bound of the variance of M at (2.1) is the variance of the linear 
regression estimator 



(/) 



M V VI = M v +dAMl -M 






( 2 . 8 ) 



where 



d l = 



[Mi 

7W, 



\{4p u (x,y)-l), 



<*2 = '^L~A ( 4 Pu(y’ z )- ] l 



fy{M 

with p n {x, y) and p n (y,Z ) being the sample analogues of the p n {x,y) and p n {v,Z ) respectively 
and f Y {M Y \ f x (M x ) and f z (M z ) can be obtained by following Silverman (1986). 

Any parametric function g(u,v) satisfying the conditions (1), (2) and (3) can generate an asymptotically 
acceptable estimator. The class of such estimators are large. The following simple functions g(u,v) give 
even estimators of the class 



g 



(1 '(w,v) = w“v p ,g (7) («,v) 



_ 1+a [u -l) 

= i - p (v - 1 y 



(3l (w, v) = 1 + a(w - 1)+ P (v - l), g ,4) (u, v) = {l -a(u - 1)- p(v - 1)} 
(u,v) = W2jM a + W,V P ,tV[ + w 2 =1 



(5) 



{6 \u,v) = au + (l-a)v p ,g (7, («,v) = exp{a(w - 1)+ p(v-l)} 



Let the seven estimators 



generated by g <l) (u,v) be denoted by My-"' = M y g^(u,v),(i = lto7). It is 

easily seen that the optimum values of the parameters a,p,W;(i-l,2) are given by the right hand sides of 

( 2 . 6 ). 

3. A WIDER CLASS OF ESTIMATORS 

The class of estimators (2.1) does not include the estimator 



M v , — M v + 



Yd 



d^Mx - M ^ ) + 1/ 2 (m z - M\\{d x ,d 2 ) 



being constants. 
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However, it is easily shown that if we consider a class of estimators wider than (2.1), defined by 



M y (g) = G 1 (m y ,u,v) 



(3.1) 



of M Y , where G() is a function of M Y , u and v such that G{M y ,l,l) = M Y and G x (M y ,l,l) = 1 . 

Gj \M y ,1,1) denoting the first partial derivative of G(-) with respect to M Y . 

/v f (7) 

Proceeding as in Section 2 it is easily seen that the bias of M Y is of the order n~' and up to this order of 

terms, the variance of M Y is given by 



Var(M Y (G) )= / [I--—. 

' 4 {f Y {M Y )Y U N ) 



i n (1 1 V 



+ 



V m n ) 



f,(M ,) 
M x f x (M x ) 



g 1 (m y ,im 



fMA 

M x f x (M x ) 



G ! (M y ,l,l)+2(4^ l (x,v)-l) 



n 


f,Gr) ( 


/,(«,) ) 


K n N y 


fz(M z )Mz l 


Mzfz(M z )} 



G, («, ,1,1)+ 2(4P n (y,z)-\)V\ 

(3.2) 

where G 2 (M Y 1,1) and G 3 (M Y 1,1) denote the first partial derivatives of u and v respectively around the point 
(M Y , (1,1). 

* 

The variance of M v is minimized for 



G 2 (M y ,l,l) = - 
G 3 (m j, ,l,l) = - 



M x f x {M x ) 
fy{M Y ) . 
M z f z (M z ) ' 
fr(M Y ) , 



(4P n (x,v)-l) 

(4P u (t,z)-1) 



(3.3) 



^ f G") 

Substitution of (3.3) in (3.2) yields the minimum variance of M y as 



min. Var 



K" 1 ) 



1 



-[ 



1 1 



4 {f Y (M Y )) 2 \m N) 
= min.Var(lVl Y (g) ) 



1 1 



1 1 



^ (4P u ( X ,y)-iy - (4P n (y, Z )-lf] 

V m nj \n N J 



(3.4) 



Thus we established the following theorem. Theorem 3.1- Up to terms of order n" 1 , 
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Var 1 



(m,' S| ) 



> 



1 1 



4{f Y (M Y f LU N) 

with equality holding if 

' f M x )m 4 



1 1 



1 1 



— ( 4 P n {y,z) — \) 



V m n _ 



\n N , 



G 2 (M Y , l,l) = -| 
G 3 (M f ,1,1) = - 



M z f z (M z ] 

fy{M Y ) 



{ 4 P n (x,y)-l) 

{ 4 P n {y,z)-l) 



J 



(G) 



If the information on second auxiliary variable z is not used, then the class of estimators M Y ; reduces to 
the class of estimators of M Y as 



M y (h) =h(m y ,u) 

where h[m Y , it) is a function of [fl Y , ll) such that H(M y ,l) = M Y and H x (M Y ,l) = 1, 



(3.5) 



8M y 



{My, 1 ) 



The estimator M Y 1 " ' is reported by Singh et al (2001). 



{ H ) 

The minimum variance of M Y to the first degree of approximation is given by 



min. Var 1 






4(/ 7 (M y )) 2 LU N) 



1 1 



V m n J 



(3.6) 



From (3.4) and (3.6) we have 



minVar(M Y ^)- min.Var(]VI Y * G *) 



n 
l n 



li I 

N Kfy(My)f 



( 4 f]i(v,r)-l)’ 



(3.7) 



/v ((j) 

which is always positive. Thus the proposed class of estimators M Y is more efficient than the estimator 
M Y H ^ considered by Singh et al (2001). 

4. ESTIMATOR BASED ON ESTIMATED OPTIMUM VALUES 



We denote 



Otj 

0, 2 



MxfxjMx) 
M Y f Y (M Y ) 
M z f z {M z ) 
Myfy{M Y ) 



( 4 P n (x,y)-l) 



( 4 P lt {y,z)-\) 



(4.1) 
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In practice the optimum values of gi(l,l)(=-ai) and g 2 (l,l )(=-a 2 ) are not known. Then we use to find out 
their sample estimates from the data at hand. Estimators of optimum value of gi(l,l) and g 2 (l,l) are given 
as 



where 



g 1 (l,l) = -d 1 

g 2 { U) = -a 2 



aj = 



mJ x ' 


K) 


Myfy\ 


k) 


M z f z \ 


^) ( 


Myf y [ 





(4^ n (x,v)-l) 



(4.2) 



(4.3) 



Now following the procedure discussed in Singh and Singh (19xx) and Srivastava and Jhajj (1983), we 
define the following class of estimators of M Y (based on estimated optimum) as 

m/ 8 ’ 1 = M Y g*(u,v,d l ,d 2 ) (4.4) 



where g*(-) is a function of (w, V,dj ,d 2 ) such that 

g *(l,l,a 1 a 2 )= 1 

dg * 0| 



g 






du 



= -a, 



(l,l,a 1 ; a 2 ) 



g 2 (l,l,aj ,a 2 ) = — — — 



g*(l,l,a,,a 2 ) = 
g*(l,l,ai,a 2 ) = 



dv 

dg * 0 



= -a^ 



(l,l 9 oci ,a 2 



dd 1 

_ dg* (•) 

5a, 



= 0 



(U,cti,a 2 ) 



= 0 



(l,l,a,,a 2 ) 



and such that it satisfies the following conditions: 

1. Whatever be the samples (S n and S m ) chosen, let M,V,d jd 2 assume values in a closed convex sub- 
space, S, of the four dimensional real space containing the point (l,l,ai,a 2 ). 

2. The function g*(u,v, dj, a 2 ) continuous in S. 

3. The first and second order partial derivatives of g * (u, V,dj ,d 2 ) exst. and are also continuous in 
S. 

Under the above conditions, it can be shown that 

e{m y {s * ] ) = M y tO^ 1 ) 
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is given by 



* (g*) 

and to the first degree of approximation, the variance of M Y 

Vat {m y [8 * ] )= min.Var(M/) 

where min.Var^My^j is given in (2.7). 

A wider class of estimators of M Y based on estimated optimum values is defined by 

My G 1 = G*(m Y ,u,v,&l jdj) 



where 



are the estimates of 



CC 2 



CX.) = 



a 2 



and G*(-) is a function of (m Y ,U,V,(X l ,d 2 j 



1) 

fyWy) 

fr[Myj 



such that 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



G * [m y ,l,l,a * ,a 2 ) = 



G* 2 {m ,a])= 

G] {m y 1,1, a*, a,) 



Gl{M Y \,\, a] , a]) 





! 

= 1 


dMy 


(m k , 1,I, a !*,«() 


0G*(-) 


* 

= — CC j 

Af y ,l,l,ct* ,cx ,2 ) 


du ( 


_0G*(-) 


* 

= _a 2 

[m y U^aj) 


dv 


_8G* (•) 


= 0 

(M r .l.l.a^.a,) 


0a* 
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G]{m Y \X^[ ,^h) 



5a* 



= 0 



I^A/y ,1,1 ,oi ! ,a 2 J 

Under these conditions it can be easily shown that 

E 



[m y (g * ] )=m y +o(/r‘) 

* (g*) 

and to the first degree of approximation, the variance of M Y is given by 

Var [m y G )= min.Var(M } <g) ) 
where min. Var [m y ° ) is given in (3.4). 



(4.9) 



(i?* I 

It is to be mentioned that a large number of estimators can be generated from the classes M Y and 

/v (g*) 

M Y based on estimated optimum values. 



5. EFFICIENCY OF THE SUGGESTED CLASS OF ESTIMATORS FOR FIXED COST 



The appropriate estimator based on on single-phase sampling without using any auxiliary variable is M Y , 
whose variance is given by 




n 

K m 



11 I 

N)4,fAM r )) 2 



(5.1) 



In case when we do not use any auxiliary character then the cost function is of the form Co-mCi, where Co 
and Cj are total cost and cost per unit of collecting information on the character Y. 

The optimum value of the variance for the fixed cost C 0 is given by 



Opt. 



Var 1 



(m,)= 



= V n 



vQ, 



l 

N 



(5.2) 



where 






4 {fr(M Y )f 



(5.3) 



When we use one auxiliary character X then the cost function is given by 

C 0 =Gm + C 2 n, (5.4) 

where C 2 is the cost per unit of collecting information on the auxiliary character Z. 

The optimum sample sizes under (5.4) for which the minimum variance of M Y is optimum, are 
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m op, = 



~oWo -riJ/Q 



n opt = 



CJVJC 2 



l(K-v l )c 1 +^. 



where Vi=V 0 (4Pn(x,y)-l) 2 . 

Putting these optimum values of m and n in the minimum variance expression of Af y in (3.6), we get 

the optimum mi n . V 3x[m y ' H ' j as 



Opt. min.Var^Mj,^') = 



k-kK+Jkcj K 



Similarly, when we use an additional character Z then the cost function is given by 



C 0 =C l m + (C 2 +C 3 )t; 



where C 3 is the cost per unit of collecting information on character Z. 

It is assumed that C 1 >C 2 >C 3 . The optimum values of m and n for fixed cost C 0 which minimizes the 
minimum variance of M y K '{pxM y° ] j (2.7) (or (3.4)) are given by 



m opt = 



' 0 \ v 0 



'o-Wl +Vi C 2 +C 3X*1 ~ V 2, 



n opt = ■ 



c„J(r,-F 2 )/c 2 +c, 



'» -r.X-, + Vic, + c, - r 2 , 



where V 2 =Vo(4Pn(y,z)-l) 2 . 

The optimum variance of M [oxM y°^ corresponding to optimal two-phase sampling strategy is 



Imin.Var 1 



(M y ^ g ^)or min.Var(M y (G *) = 



+ vl c 2 + c 3 Ao - v 2 , 



Assuming large N, the proposed two phase sampling strategy would be profitable over single phase 
sampling so long as 

Opt.Var(M y )] > Opt. min.Var(M y ^)or min.Var(M y (G *)J 
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. c 9 + c, 

i.e.— < 

C, 



V 


ly~ 
v 0 


V 


i 

kT 

l 


V 


A 


1 

i 



(5.12) 



When N is large, the proposed two phase sampling is more efficient than that Singh et al (2001) strategy if 

Opt min.Var^Mj.^jor min.Var(/V/ ) (V;i )J < Opt min.Var(M y (w) ) 



. C 7 + C 3 

i.e.— 1 < 



V, 



C, 



V\ V 2 



(5.13) 



6. GENERALIZED CLASS OF ESTIMATORS 
We suggest a class of estimators of M Y as 

3 = {M r (F) \M y [f) =f{m y ,u,v,w)) 



( 6 . 1 ) 



where U = M x / M' x , V = M' z / M z , W = M z / M z and the function F(-) assumes a value in a 
bounded closed convex subset Wc'll 4 , which contains the point (Vf Y , 1,1,1 )=T and is such that 
F(T)=M y =>F j(T)= 1, F (T) denoting the first order partial derivative of F(-) with respect to M Y around the 
point T=(M Y , 1,1,1)- Using a first order Taylor's series expansion around the point T, we get 

M y (F) = F(T ) + (m y = M y )f, (T) + (u - 1 )F 2 ( T ) + (v - 1)F 3 ( T ) + (w- 1 )F 4 ( T ) + O^ 1 ) 

(6.2) 

where F 2 (T), F 3 (T) and F 4 (T) denote the first order partial derivatives of f[m y ,U,V, wj with respect to u, 

v and w around the point T respectively. Under the assumption that F(T)=M Y and F j(T)=l, we have the 
following theorem. 

Theorem 6.1. Any estimator in 3 is asymptotically unbiased and normal. 

Proof: Following Kuk and Mak (1989), let P Y , P x and P z denote the proportion of Y, X and Z values 
respectively for which Y<M y , X<M x and Z<M z ; then we have 



M y -M y = 



1 



M x -M x = 



2 fy(My) 

1 



2 f x (M x ) 

K-m x = ‘ 

2 f x \M x ) 



(l - 2P y ) + 0 

(1-2FJ+0 



-V. 



-V. 



(i-2A)+o, 






(l - 2F z )+0 
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and 



M', - M . = 1 r(l-2F.)+0, 

2 



Using these expressions in (6.2), we get the required results. 

Expression (6.2) can be rewritten as 

M y (f) -My = [m y -My)+ (u - l)F 2 ( T ) + (v - 1)F 3 ( T ) + (w - 1 )F 4 ( T ) 



or 



M y (f) -M y = M Y e 0 + (e, - e 2 )F 2 ( T ) + e 4 F 3 ( T ) + e 3 F 4 ( T ) 



(6.3) 



Squaring both sides of (6.3) and then taking expectation, we get the variance of M Y to the first degree 
of approximation, as 



Var [m y [F) )= / 

' 4(/ r (M y )) 2 



"---1 
K m N j 



4 + 



'j__P 

v m n j 



A 0 ■+■ 



"---1 
K n N j 



A , 



(6.4) 



where 



4 = 



i + 



frfar) ' 
M z f z (M z \ 



F a 



\T) + 2(4P n (v,z)-\{ 



^00 



/l, = 



/r(4) 

.M x f x (M x \ 



fy{M Y ) 



+ 2{4P n {x,z)-\ 



F 2 2 (T) + 2(4P n (x,y)-\)F 2 (T) 

fy{M y ) ' 



MJ Z {M Z \ 



F 2 (T)F 4 (T) 



A , = 



/r(4) 

M z f z (M z ) 



fy(M Y ) ' 
M z fz{M z ), 



F 3 2 (T) + 2(4P u (y,z)-l)F 3 (T) 



+ 2 



fyjMy ) ' 

M z f z (M z ). 



F 3 (T)F 4 (T) 



The \f&v[M y F ' 1 j at (6.4) is minimized for 
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[(4^,(x,v)-l)-(4^,(x,z)-lX4^,(v,z)-l)] M s f x {M x ) 

[i-(4^,(^.z)-i) ! ] Mm,) 

= -a 2 (say) 

(6.5) 

(4^ , (*, z) - lX(4^ , (x,y) - 1) - (4 P, , (y, z) - l)(4ff , (x, z)- 1)] M z / z (M z ) 
[l-(4^,(x,z)-l) ! ] f r (M r ) 

= -a 2 (say) 



FtT , [(4^,(v.^)-l)-(4^,(.t.v)-lX4^,(x.z)-l)] MJ Z (M Z ) 

‘ [l-(4f>,(*.z)-l) ! ] ./,(«,) 

= -a 3 (say) 

* (f?\ 

Thus the resulting (minimum) variance of M Y is given by 



m inVar 1 



R (r ')= 



4 ( fy ( M ,.)) 1 



1 1 



\m N ) 



I j__ n 

y m n y 






l-(4^,(x,z)-l) : 



- + ( 4 Pn(x,y)-l) 



--zJ-W.Cv.z )-!) 2 

\n N ) 



= min.Var 



1 1 



1 



D- 



m n ) 4 (f Y (M r )) 2 1 - [4P n (x,z)~ l 2 



( 6 . 6 ) 

where 

D = [{4P n (y,z)-\)-{4P u (x,y)-l){4P n (x.z)-\)] (6.7) 

and min.Var(M y ^ G ') is given in (3.4) 

(/T) 

Expression (6.6) clearly indicates that the proposed class of estimators M Y is more efficient than the 

class of estimator M or [m y ^ j and hence the class of estimators My 1 ' * suggested by Singh et al 
(2001) and the estimator M Y at its optimum conditions. 

The estimator based on estimated optimum values is defined by 

p* = ^ : M y f = F*(m y ,u,v, w,aj,a 2 ,a 3 )| (6.8) 

where 
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£ _ [( 4 A 1 (*. y) - !) - ( 4 A 1 (*, Z ) ~ iX 4 ^ 1 l V > Z ) ~ ^4 fer ) 

[l-(4^n(x,z)-l) 2 ] fy [m y ) 

a _ (4^1 1 (x, z) - lX(4^n (x, v) - 1) ■ - (4gi 1 ( v, z) ~ l)(4gi t (x, z) - 1)] M z f z (m z ) 

lH4^Hf] f Y {^y) 



[(4 A ! ( v, 2 ) - 1 ) - (4 A , (x, y) - 1)(4 A 1 (*, z) - 1). 


] M z fz 


\m z ) 


l-( 4 Ai(x,z)-l) 2 


At 


Uj) 



are the sample estimates of ai, a2 and a3 given in (6.5) respectively, F*(-) is a function of 
[m Y , u,v, w, a l ,d 2 ,a 3 ) such that 



F*(T*) = M r 

8F * n| 

77 * (T*\ — — 



f 2 * <t*) = 

f 3 * (T*) = 

f 4 * <t*) = 
f 5 * <r*) = 
f 6 * (r*) = 

F 7 * ( T *) = 



8M 


Y 


8F * (•) 




8u 


j 1 * 


8F * (•) 




dv 


j 1 * 


0F * (•) 




8w 


j 1 * 


8F*(-) 




8a 1 


j 1 * 


8F * (•) 




da 2 


j 1 * 


0F * (•) 




da 2 


j 1 * 



= 1 

f* 

= -a x 
= -a 2 

= -« 3 

= 0 

= 0 
= 0 



where T* = (M Y ,l,l,l,ai,a 2 ,a 3 ) 



Under these conditions it can easily be shown that 

e(m ¥ ( f * ] ) = M y +0(77-') 



(6.9) 
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* ( fr*\ 

and to the first degree of approximation, the variance of M Y is given by 

Vslv{m y * ) = min.Var [m/ ) (6.10) 

where min.Var(Af j is given in (6.6). 

Under the cost function (5.8), the optimum values of m and n which minimizes the minimum variance of 
M y is (6.6) are given by 



cj(v„-v,-v s )/c, 

111 = 

- v, - v, )c, + V(r, - v, - V, )(c 2 + c, )] 

cj(r, -y,-r,)/c 2 

n ° P ' t M ~ v < ~ )c, + V(r, -r 1 + r, Xc 2 + c, )] 

where 



( 6 . 11 ) 



^3 = 



D V n 



\-{4P u {x, Z )-l) 2 



for large N, the optimum value of min. Var(/h/ j/ F) ) is given by 



Opt. 



min.Var 1 



{m y [f) \ 



#0 -V, -V,)c x +J{v l -v 2 + v 3 \c 2 +c 3 ) 



c n 



( 6 . 12 ) 



(6.13) 



The proposed two-phase sampling strategy would be profitable over single phase-sampling so long as 

Opt. Var (m y )] > Opt. min.Var(M Y (/ J )j 



C 2 +C, 

i.e. — -< 






It follows from (5.7) and (6.13) that 

Opt. min.Var(M K l/ ! ) < Opt. inin.Varj/V/j/ 7 ) 



if 



ck 

1 

1 

1 

1 




\C 2 +C 3 


1 Vl C 2 


l ^-V 2 +V 3 j 




J Q 


ur,-r 2 +v,)c, c, 



(6.14) 



(6.15) 



for large N. 

Further we note from (5.11) and (6.13) that 
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Opt. 



min.Var^Mj ,^|]< Opt. min.Var(M y ^^orM y G )J 



if 



c 2 +c 3 

c i 



VCo-r,) 




K- 


Vi-V>) 


_ yi 




+ ^3 )“V 


' y , - v 2 



(6.16) 
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1. INTRODUCTION 

The discrepancies between the values exactly obtained on the variables under consideration for sampled 
units and the corresponding true values are termed as measurement errors. In general, standard theory of 
survey sampling assumes that data collected through surveys are often assumed to be free of measurement 
or response errors. In reality such a supposition does not hold true and the data may be contaminated with 
measurement errors due to various reasons; see, e.g., Cochran (1963) and Sukhatme el al (1984). 

One of the major sources of measurement errors in survey is the nature of variables. This may happen in 
case of qualitative variables. Simple examples of such variables are intelligence, preference, specific 
abilities, utility, aggressiveness, tastes, etc. In many sample surveys it is recognized that errors of 
measurement can also arise from the person being interviewed, from the interviewer, from the supervisor or 
leader of a team of interviewers, and from the processor who transmits the information from the recorded 
interview on to the punched cards or tapes that will be analyzed, for instance, see Cochran (1968). Another 
source of measurement error is when the variable is conceptually well defined but observations can be 
obtained on some closely related substitutes termed as proxies or surrogates. Such a situation is 
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encountered when one needs to measure the economic status or the level of education of individuals, see 
Salabh (1997) and Sud and Srivastava (2000). In presence of measurement errors, inferences may be 
misleading, see Biemer et al (1991), Fuller (1995) and Manisha and Singh (2001). 

There is today a great deal of research on measurement errors in surveys. An attempt has been made to 
study the impact of measurement errors on a family of estimators of population mean using multiauxiliary 
information. 

2. THE SUGGESTED FAMILY OF ESTIMATORS 

Let Y be the study variate and its population mean u 0 to be estimated using information on p(>l) auxiliary 
variates X b X 2 , ...,X p . Further, let the population mean row vector Jut ' = (|U , , \X 2 , • • •, \X p ) of the vector 

X ' = (x t , X 2 , X p ] . Assume that a simple random sample of size n is drawn from a population, on the 

study character Y and auxiliary characters X b X 2 , ...,X p . For the sake of simplicity we assume that the 
population is infinite. The recorded fallible measurements are given by 



yj= Y j +E j 

x o = x ij +, V i = — 

j = 1,2, •••,«. 

where Yj and Xy are correct values of the characteristics Y and X ; (i=l,2,.„, p; j=l,2,..„ n). 

For the sake of simplicity in exposition, we assume that the error Ej's are stochastic with mean 'zero' and 
variance a (0) 2 and uncorrelated with Yj's. The errors Py in xy are distributed independently of each other 
and of the Xy with mean 'zero' and variance cpi) 2 (i=l,2,...,p). Also E/s and py's are uncorrelated although 
Yj's and Xy's are correlated. 

Define 
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U, = —,(/' = 1 , 2 , •••,/?) 



T 

U = I 



= (u 1 ,u 2 ,—u p ) Up ,e T =(l,l,-",l) 1X p 



1 /? 

v =-!>,. 
n j= \ ■ 

i n 

*i =~ll x u 
n j= i 



With this background we suggest a family of estimators of ju. 0 as 

i i g =g(yy) 



(2.1) 



where g\V,U 



(v,w r ) 



is a function of y, u l ,u 2 ,---,u 



such that 



g (u 0 ,e T ) ^ 



8y 



= 1 



^o. e 



and such that it satisfies the following conditions: 

1. The function g(v,U T ^j is continuous and bounded in Q. 

2. The first and second order partial derivatives of the function g(y, U 7 j exist and are continuous and 
bounded in Q. 

To obtain the mean squared error of |i , we expand the function g{y, U 1 ) about the point (lio C 1 ) in a 



second order Taylor's series. We get 



= 



= g(flo^ r )+(v-fl 0 } 






dy 



+ 



{u-e) T g (l \^,e r ) 



{Vo> e 






dy 2 



+ 2{y- V-oh-e) 



8y I 



+ (u - ef g [2) {y*,u * T \u~ e)) 



(2.2) 
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where 



y*= Ho + © ( v — J u 0 ), u * = e + 0(u-e),0<Q <1; g (1, (-) 

denote the p element column vector of first partial derivatives of g(-) and g (2) (-) denotes a pxp matrix of 
second partial derivatives of g(-) with respect to u. 

Noting that g(p 0 ,e T )= Po, it can be shown that 

Po +0(n x ) 



(2.3) 

which follows that the bias of fl is of the order of n" 1 , and hence its contribution to the mean squared 



error of ju, g will be of the order of n' 2 . 

From (2.2), we have to terms of order n , 

MSE (^g ) = E \y - Ho ) ■ + ( U - e Y g ( Xo,e T ) f 

= e[{v = p 0 ) 2 + 2 ( V - Ho ){u - e) T g ( \y) 

+ (g ( % 0 y)) r (u-e)(u-e) r (g ( % 0 y)} 

= MvoiCo 2 +C ( l))+2v 0 b r g%S)+(g% 0 stf 

i) 



where b T =(bi,b 2 ,...,b p ), h^poAC^^l^, ...,p); 
C,=a,/p„ C (i) = a/p,, (i=l,2, ...,p) and C 0 =c 0 /po, 





c 2 +c 2 } 


p l2 C l C 2 


Pl 3 QC3 


•• Pi P C,C P 




Pi 2 Ci c 2 


C 2 +C( 2 2) 


p 23 C 2 C 3 


•• P 2p C 2 C p 


A = 


p 13 C[ c 3 


p 23 C 2 C 3 


c 3 2 + c ( 2 ) • 


•• Ps p C 3 C p 




_P lp C l C p 


P 2 p c 2 c p 


P 3p C 3 C p ■ 


■■ c ; +c l ) 



The MSE(|rt g ) at (2.4) is minimized for 

g ( ‘W) =-[i 0 A~ l b 



(2.4) 



(2.5) 
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Thus the resulting minimum MSE of JLX „ is given by 



min.MSE(|lJ= (n 0 2 /n)[C 0 2 +C 2 0) -b T A x b\ 



( 2 . 6 ) 



Now we have established the following theorem. 
Theorem 2.1 = Up to terms of order n" 1 , 



MSE(H a )> (j , 0 2 !n{cl +C ( ! 0| -b T A-'b\ 



(2.7) 



with equality holding if 



g [ \ 0 ,e T )=-\i 0 A~ l b 



It is to be mentioned that the family of estimators |ll at (2.1) is very large. The following estimators: 



I i, w = vl>, E®, = 1. 



i= 1 






J 



i = 1 



[Olkin (1958)] 



, . p f % "N p 

lV — 7 Z 03 ' -1 ’ [Singh (1967)] 



i= 1 



Vi) 



i = 1 



( 3 ) _ TT M 






= y 






Z 03 / = !» [Shukla (1966) and John (1969)] 



1=1 



1=1 



( 4 ) - J=1 



Z 03 ,^, 



Ee = V 



Z 03 ^ 

1=1 



Z* 33 / = [Sahai et a/ (1980)] 



^„ (5) = vn ~ 7 Z°T = 1, [Mohanty and Pattanaik (1984)] 






1=1 



\ x i J t = t 



~ ( 6 ) - 

^ = y 



fpavY 1 JL 



Z — '~~ L , Z* 33 = 1’ [Mohanty and Pattanaik (1984)] 

/=! H/ J /=! 
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^’=?n 



1=1 



r ~ \® 

^i_ 



y 

Xco ; = 1, [Tuteja and Bahl (1991)] 



i=i 



-|-1 



- (8) - 
P, = V 



y 

I 



<"+ 



Xco,. =1, [Tuteja and Bahl (1991)] 



;=i 



- ( 9 ) 

^ =>' 



- (10) - 
= v' 



°vi +Z®/ 

1= 1 



°vi +Z®< 

i = 1 





\ 


i+ 




v + 


)_ 


y - 


w 




- 



M 



p+i 

Ecoi =!• 

i=l 



P+1 

» Z®«- =1 - 

i= 1 



- (11) - 
f- 1 , = v 



s 

i=l 



CO, 



V 

V*</ 



p 

+ x 






<= 4+1 



r * p v 1 

Z®/ + Z ro 

i=i i= ? +i y 



; [Srivastava (1965) and Rao 



and Mudhalkar (1967)] 



K v II (a,, 's are suitably constants) [Srivastava (1967)] 






1=1 



I h) 






^ - 77 



?n 



i=i 



' 


/ - \ a < 1 




[ +_ 




2- 











^ [Sahai and Rey (1980)] 



|a 



0 4 ) _ 77 



>n 



{(a, +a. i (x i - |aj} 
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etc. may be identified as particular members of the suggested family of estimators JLl . The MSE of these 
estimators can be obtained from (2.4). 



It is well known that 

v(y)=(^/»Xc„ 2 + cj,) 

(2.8) 

It follows from (2.6) and (2.8) that the minimum variance of |.l ( , is no longer than conventional unbiased 
estimator y . 

On substituting ct ( 0 ) 2 =0, cj (i) 2 =0 Vi=l,2,...,p in the equation (2.4), we obtain the no-measurement error case. 
In that case, the MSE of JLl , is given by 

MSE(jI g )= i[c 0 > 0 2 +2 MV (1 W) +fe* (I W)) r ^*t * W (m.^)). 

= MSE(jI g *) 



(2.9) 

where 






= g* 



... ^ 
v E, E 2 E 



p J 



(r.un 



(2.10) 

and Y and X ; (/ = 1,2, •••,/;) are the sample means of the characteristics Y and X; based on true 
measurements. (Yj,X;j, i=l,2,...,p; j=l,2,. . .,n). The family of estimators |4 g * at (2.10) is a generalized 
version of Srivastava (1971, 80). 

The MSE of jl „ * is minimized for 



g 




-A *- 1 bn 0 



( 2 . 11 ) 

Thus the resulting minimum MSE of ji „ * is given by 
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min.MSE^ *) = — [c o 2 - b T A * _1 b] 
n 



2 




(2.12) 

where A*=[a*jj] be a pxp matrix with a*;j = pyCiCj and R stands for the multiple correlation coefficient of 
Y onX b X 2 ,...,X p . 

From (2.6) and (2.12) the increase in minimum MSE (jl,, ) due to measurement errors is 



obtained as 



/ 2 \ 

min.MSE^J-min.MSE^ *)= [c^+b 7 A* ' b -b T A l b\ 



V 

>0 



This is due to the fact that the measurement errors introduce the variances fallible measurements of study 
variate Y and auxiliary variates X;. Hence there is a need to take the contribution of measurement errors 
into account. 



3. BIASES AND MEAN SQUARE ERRORS OF SOME PARTICULAR ESTIMATORS IN THE 
PRESENCE OF MEASUREMENT ERRORS. 



To obtain the bias of the estimator |i , we further assume that the third partial derivatives of 
also exist and are continuous and bounded. Then expanding g(y, U ' ) about the point 
(v,ll 1 ) = (|U, 0 , e T ) in a third-order Taylor's series we obtain 

\%0I 



v* = S 



(Ho,e r )+(v-H 0 )- 



dy 



+ (u-e) T g ( 1 \v 0 ,e r ) 







1 

+ — 

2 



fc v-^y 8 -W 



dy z 



+ 



^{y - v Q X u ~ e Y g { \^ T ) 






+ {u-e) T (g { '% 0 ,e T )Xu-e)) 
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(3.2) 

where b T =(bi,b 2 ,. . .,b p ) with bi=p O iC 0 Ci; (i=l,2, , p). Thus we see that the bias of ji depends also upon 

the second order partial derivatives of the function on g[v,U T ) at the point (po,e T ), and hence will be 
different for different optimum estimators of the family. 

The biases and mean square errors of the estimators i = 1 to 1 8 up to terms of order n" 1 along with 

the values of g (1) (|i 0 ,e T ), g (2) (p 0 ,e T ) and g <12) (p 0 ,e T ) are given in the Table 3.1. 
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Table 3.1 Biases and mean squared errors of various estimators of (i 0 
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Table 3.1 Biases and mean squared errors of various estimators of (i 0 







4. ESTIMATORS BASED ON ESTIMATED OPTIMUM 



It may be noted that the minimum MSE (2.6) is obtained only when the optimum values of constants 
involved in the estimator, which are functions of the unknown population parameters p 0 , b and A, are 
known quite accurately. 

To use such estimators in practice, one has to use some guessed values of the parameters p 0 , b and A, either 
through past experience or through a pilot sample survey. Das and Tripathi (1978, sec. 3) have illustrated 
that even if the values of the parameters used in the estimator are not exactly equal to their optimum values 
as given by (2.5) but are close enough, the resulting estimator will be better than the conventional unbiased 
estimator y . For further discussion on this issue, the reader is referred to Murthy (1967), Reddy (1973), 
Srivenkataramana and Tracy (1984) and Sahai and Sahai (1985). 

On the other hand if the experimenter is unable to guess the values of population parameters due to lack of 
experience, it is advisable to replace the unknown population parameters by their consistent estimators. Let 

(f) be a consistent estimator of (j)=A -1 b. We then replace 4> by (f) and also p 0 by V if necessary, in the 

optimum ju g resulting in the estimator £l „( est ) , say, which will now be a function of V , u and <|>. Thus we 

define a family of estimators (based on estimated optimum values) of p 0 as 



M g (es,) =g**(y,U T ,0 T ) 



where g**(-) is a function of 



[y,u T ,<j> T ) 



such that 



g**(/u 0 ,e T ,</> T )= /u 0 for all ju 0 , 
dg* *(•)! 



dy 



= 1 



[n 0 ,e T ,/) 



(4.1) 



dg * *01 



_ 3g0 



du 



U.«V) du 



= ~M 0 ^ l b = -/uj 



(uo, e ) T 



(4.2) 




and 



8g * *(■) 

d(j> 



= 0 



With these conditions and following Srivastava and Jhajj (1983), it can be shown to the first degree of 
approximation that 

MSE (Ar(«o) = min -MSE(// g ) 



f 

V n J 



po + E’(o ) -b T A~ x b\ 



Thus if the optimum values of constants involved in the estimator are replaced by their consistent 
estimators and conditions (4.2) hold true, the resulting estimator jU g ( est ) will have the same asymptotic 



mean square error, as that of optimum jU a . Our work needs to be extended and future research will explore 
the computational aspects of the proposed algorithm. 
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The purpose of this book is to postulate some theories and test them 
numerically. Estimation is often a difficult task and it has wide application in 
social sciences and financial market, in order to obtain the optimum efficiency for 
some classes of estimators, we have devoted this book into three specialized 
sections. 
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